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Abstract 



We show that a Galilean invariant version of fluid dynamics can be derived by the methods 
of statistical dynamics using Maxwell's balance equations. The basic equation is non-local, and 
might replace the Boltzmann equation if the latter turns out not to have global smooth solutions 
■ in general. As an approximation, a local form of the equations of motion is derived. It turns out 

to be a version of the Navier-Stokes system, obeying the Stokes relation, and with the viscosity 
coefficient rising as O 1 / 2 with temperature 0. The new feature is the presence of the Dufour 
effect for a gas of a single component. This ensures that the principal symbol of the parabolic 



system is non-singular. 

o 



1 Introduction 



A central problem for mathematical fluid dynamics is the derivation of the Navier-Stokes equa- 
tions (N-S) starting from a reversible dynamical theory such as classical hard spheres or quantum 
mechanics. To achieve this, it seems essential to know that the N-S equations themselves possess 
smooth solutions for all times, for a large enough class of smooth initial values for the fields. Hence 
the latter problem is revealed as the key question, for one version of which a Clay Millenium prizes 
is offered. The prize version describes an incompressible liquid under isothermal conditions. Nash 



pi| , 32] had shown that, given smooth initial conditions, there exists a unique smooth solution for a 



small enough time. The question of smooth global-in-time solutions remains open except for small 
initial conditions. To model driven systems such as Benard convection, the isothermal condition 



must be relaxed; recent numerical studies [2C] show qualitative agreement with experiment. In [20], 
the condition of incompressiblitiy, divw=0 is maintained, and the energy equation is modified by 
the addition of a bouyancy condition. This is expressed by requiring that the liquid in hot regions 
is less dense than in cold regions, a rather ad hoc procedure. We shall argue in Sect. (2) that N-S is 
at the boundary of a more regular class of models. In N-S, the pressure is infinite, but its place is 
taken by a surrogate pressure determined by the requirement of self-consistency. Thus in N-S, the 
pressure is a balancing item like petty cash, much beloved by accountants, which can be adjusted 
to cancel errors made elsewhere in the calculation. 

In this paper, C-N-S- T will denote the system of five coupled non-linear partial differential 
equations known |^7j as 'compressible Navier-Stokes with temperature'. This might or might not 
be an easier problem than N-S, but it is certainly more widely applicable. A difficulty with C-N- 
S-T is that the symbol of the elliptic operator is singular. To derive N-S, some authors divide the 
problem into two parts ||; first, to show that reversible dynamics is well represented in some limit 
by a stochastic process; then to show that the Fokker-Planck equations of this process gives rise 
to the N-S equations. The latter is only partially achieved in Q. Other authors start with some 
version of stochastic dynamics [^, J5J, and prove things; this approach will be adopted here, 
as it avoids the much harder first part. We shall adapt information geometry ^ |25| to the 



dynamics of a rare gas. Thus, the state of the system is not a point in phase space, but a measure 
H on it, and the dynamics is a path in the space E of measures. In the usual treatment j|] the 
currents of the conserved variables are given exactly by expressions involving higher moments of 
the same random fields; the time derivatives of these higher moments involve yet higher moments. 
The whole system goes on for ever, and is called the BB GKY hierarchy. Inasmuch as the system is 
equivalent to classical mechanics, it is reversible and shows no dissipation. Artfulness is needed to 
'close' the system in terms of the slow fields; that is, to write the time-derivatives of the slow fields 
in terms of themselves, by truncating the system. It has proved possible to get a variety of kinetic 
equations from the hierarchy, by taking a limit appropriate to the physical situation ||. These 
limiting systems exhibit dissipation. It then remains to show that the solutions to the limiting 
system are limits of solutions to the BB GKY hierarchy. This programme has proved to be difficult 
to complete. 

Information dynamics offers an alternative. In the case studied here, where the potential be- 
tween the particles is zero outside the hard core, the states in local thermodynamic equilibrium, 
LTE can be computed. The LTE states are products over the lattice; the state at a site x is of the 
form 

K^x) = N x p(x,k), 

where p is Maxwellian. The set of such states make up the information manifold, A4. Any state 
/i£E, having finite means for the slow variables, has a reduced description, denoted nQ\ here Q 
is the the non-linear projection onto M, acting on the right, which maps \x to the state in A4 with 
the same means for the slow variables as \x. The Gibbs principle 23, ^] states that \iQ is the state 
of maximum entropy having these means. In the original formulation of information dynamics [21], 
in the time interval (0, t) the state \x evolves under the reversible dynamics of classical mechanics, 
to /i(t) say. This state has a much simpler description by the LTE state [i(t)Q, which cannot be 
distinguished from fj,(t) by measuring the slow variables. The orbit {fj,(t)Q : t > 0} in M. was 
intended as the thermodynamic evolution. It is clear that the entropy of fJ,(t)Q is not less than that 
of there is a transfer of information into inaccessible degrees of freedom by the reversible motion 
[]|]. It is not always true that entropy increases along the orbit, as is seen if the classical motion were 
periodic. Another version of information dynamics was adopted in Q; there, the reversible motion 
took place for a very small time t, and the reduced description fi(t)Q was used, instead of fj.(t), 
as the initial state of the next step. This gives a discrete-time semigroup, with increasing entropy; 
however the time-step cannot be taken to zero, without sending the rate of entropy production to 



zero as well [25, |4|. It is necessary to keep the time-step positive; it represents the relaxation time, 
and the map Q implements the thermalisation of the state n(t). The challenge is to do this in a 
way that is invariant under the Galilean group (denoted by Q below). 

Information dynamics has been extended [^] to allow stochastic dynamics; then the time-step 
can go to zero, still giving a non-zero rate of entropy production. Another idea is to allow state- 
dependent transition rates [Q]. With these changes, one may call the theory statistical dynamics. It 
is designed to obey both the first and the second laws of thermodynamics, but otherwise puts few 
constraints on the form of the dynamical equations. The choice made for the dynamics determines 
the nature of the system under discussion. At first sight, statistical dynamics has too much noise; 
a simple application is shown 37] to lead to mass diffusion and the Soret effect for an inert gas at 
rest, contrary to the literature ||. Indeed, without a velocity field, the theory 'has not got off the 
ground' [26|. Truesdell [38] ironically says "results of this kind are described by kinetic theorists as 
'corrections to hydrodynamics' ". 

In this paper, we apply statistical dynamics to the case of an inert gas of a single type. We 
arrive at C-N-S-T, but with one extra term, a Dufour effect. Thus the intuitively attractive 'method 
of Maxwell', [29| in which we compute the gain and loss of particles in a small time interval at 
each point x, is successful. The new idea is to postulate that in the state some but not all the 



particles are thermalised; those that are, are described by the LTE state JL. However, this state is 
NOT equal to [iQ\ By construction, our model is the the continuum limit of a non-linear Markov 
process on a lattice, with a bistochastic transition matrix. This might be a possible starting point 
for a proof that smooth solutions exist. The method of Maxwell is treated in 14], p 93, and in 
Balian ||, but it is abandoned as too primitive, in favour of the Boltzmann equation. This might 
have been too hasty, especially if the BE turns out not to have smooth global solutions. 

In Sect. (2) we study a discrete model of hard spheres, for which the equilibrium state exactly 
factorises. Sect. (3) contains a discussion of the idea that the true state \x should be distinguished 
from its thermalised part, /Z, using the analogy of aircraft in an airfreight company. The hopping 
rules of the nonlinear Markov chain are presented, and related to the mean free time. We also 
suggest a (/-invariant collision function. From this, we get an explicit formula (|4(]) for the mean 
free time, t£. In Sect (4) we start with the fundamental relation (^4|) expressing the full state as 
a non-local integral of the thermalised state. The dynamics is expressed in terms of the BBGKY 
hierarchy. We obtain the Euler equations as the zero th approximation, and obtain a useful short 
version of these. In Sect. (5) we find the differences of the means of the slow variables in the 
states /i and Ji, and show how the method can be extended to other variables. This enables us to 
compute the BBGKY moments in terms of the means in the full state and thus arrive at a version 
of C-N-S-T. This exhibits the Dufour effect, contrary to the literature M. 



2 The Thermostatics of an Inert Gas 

We take space to be A C (aZ) 3 , and suppose the length a, representing the diameter of a molecule, 
to be so small compared with the variation of the macroscopic fields that we can replace all sums 
over A by integrals. The possible configurations of the fluid are the points in the product sample 
space 

xeA 

so a configuration is specified by the collection {tox}xeA- For each x, 



^ = {0,(eZ) 3 } 



Here, e is a small parameter having the dimension of momentum; for example, we could take oe = h, 
the semi-classical division of the phase-space of a particle into cubes of volume h 3 . If the system 
is in a configuration uj, such that ujx = 0, then we say that the site x is empty. If lux = k, we 
say that the site x is occupied, by a particle of momentum k. This simple exclusion of more than 
one particle on each site incorporates the hard-core repulsion between the particles, which are thus 
hard spheres sitting at some of the points of A. The field point of view enables us to avoid the 
Gibbs paradox. 

The state of the system is a probability on f2, denoted by fi. We denote the set of states by S. 
The 'slow variables' of our model are the 5 extensive conserved random fields 

\r ( \ f Mu x = ® 

N *M = \ l Mu x = k (1) 

f < \ J if ^ = , . 

Cx{UJ > ~ y k • k/2m + <S>(x) if u x = k {1 

i s J Moj x = 
\ k it uj x = k 



Here, &(x) is the external potential energy per particle. The variables appearing in the C-N-S-T 
equations are simply related to the mean fields in the state \i: 

N x = E M [jVx]; E x = -E tl [E x ]; U X = E^[V X }. (4) 

In information geometry, the specification of the slow variables determines the information man- 
ifold M, which in the context of fluid dynamics consists of states in LTE (local thermodynamic 
equilibrium). Such a state is specified by five canonical fields, dual to the mean fields: [3 X , (s, Cxi 
and has the form 

H{u) = J] Ex 1 exp {-ZxN x (u) - (3 x 8 x {u) - £ x • V x (lu)} . (5) 
xeA 

In finding the partition function 



/2irm\ 3 / 2 

1 + £_ \~Px~) e ^{-tx-Px<5>{x)+mC x -C x /Wx} (6) 



we have replaced the sum over the momentum lattice of size e by a Gaussian integral. The product 
structure of an LTE state means that an observable at a point of A is independent of an observable 
at any other. The state fi can be written in Maxwell form 

fi = N xP (x, k) = N X Z- 1 exp {-(3 x ${x) - (3 x k • fc/(2m) - C* ■ k} , (7) 

where 



_ 3 /27rm\ 3 / 2 f mC x -C x \ 



(8) 



We note the identity for each x 

3 = 1 + e" ? Z. 



The external potential does not influence the local velocity distribution, as it is cancelled out by 
the partition function. The mean fields (||) are related to the canonical fields by 

Ex = ~lv^x = N{x)Ux) + ^ + ^%^] (9) 



d/3 x & v ' V w Wx 2(3, 



>x 



d Hjb - 1 Ze^a? , % 

iVa; = -^ — ogHa; = — = (10) 

U x = --lo g ,x = (11) 

The formalism breaks down if (5 is zero or infinity, or if N vanishes, but the case of a fluid at rest, 
II = 0, is within the information manifold, M. 

Historically, the intensive variables used in the N-S equations were the chemical potential —£//?, 
the velocity field u = — £/(3 and the temperature = (k B (3)~ l . We shall eliminate £ in favour of 
the mass-density p = a~ 3 mN using fllO|), which leads to 



e-te =Z X 1 N X /(1-N X ). (12) 
The mean occupation per site N x obeys < N x < 1. The (von Neumann) entropy of any state \x 

S(M):=-* B X>(w)logM(w)- (13) 



is 



Gibbs knew that the state of maximum entropy, among all states with the given means of the total 
energy and number of particles, is the exponential state that he called the grand canonical state 
(23, 22, p5j ]. This is a simple exercise in Lagrange multipliers. If the mean fields depend on x, then 
the state of maximum entropy has the same form, in which the canonical fields £, (3 and £ now 
depend on x. 

In this section we study the system in equilibrium, and denote by E, N and II the total values 
of the energy, number and momentum; then (^) gives for the entropy 

OS(fi) = E + k B G£N -u-Il + k B eiogZ. (14) 

Compare this with the thermostatic formula 

OS = E + k B 0£N -u-U + PV (15) 

(note that the term u • II is omitted in [27|, eq. (1.17)), where P is the pressure and V is the 
volume; we see that 

P = k B log 3 = k B 9a- 3 log H. (16) 

If there are N = J2x Nx particles, and Vo is the smallest volume they can occupy (one per site), 
then V = a 3 N and N x = V /V. Also, 

Ha; = (1 - Nx)' 1 = 1 + V /(V - Vq). 

Thus at equilibrium, we have the equation of state 

„ *.e„„. J, , v« 



- + < 17 > 

For small Vq/V this is close to the van der Waals gas 

( P + A /V 2 ){V-Vo) = Nk B S (18) 
with ^4 = 0. Unlike the case A > 0, this model shows no failure in convexity in its isothermals. 



3 The Statistical Dynamics of the Gas 
3.1 An Airfreight Model 

Consider a gas of free particles in a box with reflecting walls, in equilibrium; then u, and p do not 
depend on x. There is still a lot going on. In a volume d 3 x around x, a particle of momentum k, 
which is present with probability Np(k), moves in the direction of the unit vector k, to be replaced 
in time t by a particle with the same k arriving from the point x — kt/m. This replacement 
was present with exactly the same probability. The larger k is, the further away is the source of 
the replacement. In this picture, equilibrium is described by a huge game of musical chairs; only 
the indistinguishability of the particles prevents this from being detected. Now look at the same 
mechanism, but where p, II and O depend on x. That is, we now consider the Knudsen gas. There 
is no longer exact replacement of the lost particles at x; the parameters p . . . change with time. It 
might seem that the system gets closer to equilibrium, since the parameters start to become more 
and more nearly constant. Of course, the entropy is constant in time, as the system is Hamiltonian, 
(free, even). The apparent increase in entropy associated with the slow variables is exactly matched 
by a reduction of entropy in inaccessible observables ||. Thus, taking the initial state to be in 
LTE, the random variables M ... had independent values at every point. But after some time, 



the free motion introduces correlations between very far points; if x at time t has a particle with 
momentum k, then x — kt/m must have had a particle of momentum k at time 0; so it did not 
have a particle of momentum kl ^ k. Thus x + (k' — k)t/m has no particle of momentum k' at time 
t, a statement correlated with the assumption, above, about x at time t. Correlations like these 
might at any time show up in behaviour quite unlike that of a system near equilibrium. For spin 
systems, the spin-echo effect is such an example [|. In a real gas, we do not expect any surprises 
such as the spin-echo effect. This is due to the interactions, idealised by collisions, which remove 
correlations between distant points, and also help to redistribute energy and momentum. Collisions 
do not contribute to transport; on the contrary, they inhibit the free flow: the diffusion constant is 
inversely proportional to the collision cross-section. 

The concept of whether a system is thermalised or not is independent of the Galilean frame 
of reference used in the description. This is expressed mathematically by the fact that the set of 
equilibrium states is mapped to itself by the group Q. The set M. of states in LTE is also mapped 
to itself by Q, which we interpret as saying that the concept of partially thermalised systems is also 
invariant under Q. Physically, a gas consists of some (most) particles that are thermalised, and 
are described by an LTE state with means equal to the averages over the thermalised particles. 
A smaller number are not well described in this way; in particular, their correlations with other 
particles are underestimated by assuming LTE. Moreover, a particle that is thermalised at time t 
will move under its free motion to regions at different density and temperature, and so after some 
time (how short depends on the gradients) it will not be well described as being thermalised. On 
the other hand, particles left out of the count of thermalised particles make collisions during their 
relaxation time, and return to the thermal state. The mass, energy and momentum will on average 
be conserved, but a particle can leave the LTE state at one point, and another can return at another 
point, so the microscopic currents describing all particles may differ slightly from the currents of 
thermalised particles as described by the LTE state. This division of the state into thermalised 
particles and the rest differs from the division presented e.g. in H, p. 160, where the LTE part of 
the state gives the same expectation values as the true state. On the contrary, in our division, the 
state JL does not give the same means as the true state fi. A different division can be found in |18| , 
p 229, where nF[ is taken to be 'the density of particles that are not undergoing collisions at the 
given instant'. From now on, the variables N,E,H, etc., refer to the full state, and written with 
bars, they refer to the LTE state of the thermalised particles. 

The dynamics of a classical gas of hard spheres is similar to that of an airfreight company, 
whose planes fly between airstrips arranged in a lattice A. In calculating the overall transport of 
goods, the company uses statistical methods; they have records only of the local averages, at each 
airstrip, of the number of planes, their velocities, and their kinetic energies, at time t = 0. Every 
plane is instructed to select a velocity from the Maxwell distribution at its airstrip, and to fly with 
this velocity until it meets another plane in its airspace. Both planes must then land very briefly, 
and record their presence, momentum and energy to the local computer. This recalculates the 
updated values of N, II, E for this strip, and instructs them to take off with new velocities drawn 
from the updated Maxwellian. The problem is to find a theory which can predict the average 
transport of goods, M, V, £, without reading the local computers. To account for the transport 
of N, II, E, we must introduce accounting system non-local in space and time: we know that just 
after a landing and take-off, the distribution of velocities is Maxwellian. This simple fact leads us 
to the fundamental equation (^4|). 

We want the dynamics to satisfy the second law of thermodynamics. This is ensured in a 
Markov chain if the transition matrix is bistochastic. Physically, if the gas has no velocity u, the 
transition rate from x to y by a particle of momentum k is the same as the rate from y to x by a 
particle of momentum —k. This expresses time-reversal invariance of the transition matrix. That 
is, if 7" : lj i — ► cut is the time-reversal map on the sample space, we say a transition matrix T obeys 



time reversal symmetry if 



T(oj\u}') = T(u>'t\lut) for all uj,uj' £ fl. (19) 

We note the following lemma, whose proof is simple: 

Lemma 1 Let T be a stochastic map obeying then T is bistochastic. 



Our model dynamics will be given by a stochastic map obeying (|19|), and so, by the lemma, will 
be bistochastic, and so entropy-increasing. 

3.2 The Hopping Rules 

In this section, we give hopping rules for the case of zero external field, <I> = 0. The dynamics will 
be invariant under Q. We start with a model in discrete space-time. The discrete dynamics will 
be given by specifying a hopping probability in one time-step. In the classical hard-sphere model, 
it is to be expected that on average, particles of different speed travel through the same amount 
of material before thermalising. Let £(x, k) denote the average distance travelled by a particle 
starting at x with momentum k; to begin with, assume that k lies along one of the basis vectors 
of the lattice. £ is called the mean free path, and it generalises an idea going back to Clausius [16|; 



it is going to be the mean of a random distance r, the free path between collisions. We assume, 
as part of the model, that the particle thermalises on its first collision. The relaxation time t of a 
particular particle depends on its speed; for a particle travelling the free path r, t is the time taken, 
rm/\k\. We therefore cannot choose a unique time- interval for the time step of all processes, and 
it seems difficult to construct a discrete-time stochastic process. We overcome this complication 
by noting the rate \k\/(rm) at which the process transfers mass, energy and momentum; we can 
then move to a continuous time process with the same rate. We shall work with t and its mean, 
te, rather than with the free path r and its mean, I; tg has the advantage of being the same in all 
inertial frames. 

The dynamics must be such as to conserve the totals 

M ■■= E Mxi £ ■= E £ x; v := E 

xeA cceA xeA 

These random variables divide f2 into simultaneous level sets, the mass-shells, energy-shells, and 
momentum shells, thus: 

n N,E,U : = 6 Q : -^M = N > £ ( UJ ) = E i V ( UJ ) = n | • 

Clearly, 

n = U n N,E,n- 

The dynamics, the Markov matrix T, must be chosen so that a point uj jumps to another point in 
the same shell. We cannot expect this to be given by a symmetric Markov transition matrix: the 
inverse process involves a change of sign for k; however, we shall be able to construct a suitable 
bistochastic map. In fact, our Markov transition matrix T will be a convex mixture of permutations 
that move a configuration u\ to U2, where in uj\ there is a particle at x and a hole at y, and in LU2 
the opposite holds, with the mass, energy and momentum that was at x transported to y. This 
move is only possible if all the points between x and y are empty. Moreover, to ensure that at 
the end of the flight the particle returns to the thermalised fold, the site one place past y must 
be occupied. We postulate that a particle moves in a straight line along empty sites until it meets 
a filled site; it then thermalises at the last empty site, y say, and dumps its mass, energy and 



momentum there, which joins the mass, energy and momentum of the state \xy. We postulate one 
such jump for each r = sa, where s is an integer, and for each momentum k at x, and then for 
each x. 

Suppose that the site x is occupied, with momentum pointing in the direction of one of the basis 
vectors of the lattice, say k = \k\e. The probability that the sites x + s'ae be empty, 1 < s' < s, 
and the site x + (s + l)ae occupied, is 

s 

1] (1 - N x+S > a e) N x+{s+ i)ae (20) 

s'=0 

One can check that e.g. if Ny > is independent of y for large enough \x — y\, then 

s 

E II C 1 - N x+s'ae) N x+{s+1)ae = 1. (21) 

s s'=0 

This just expresses that with probability one, the number of holes adjacent to x on the line joining 
x to infinity along the direction e must be some integer. The product in ( pOj ) is a marginal 
probability of the state /x, and so is linear in /i; however, it is a polynomial of indefinite degree in 
the variables N x , in terms of which the equations of motion are to be written. In discussing the 
flow of mass, energy and momentum at the point x, it is convenient to include this factor in the 
hopping probability, rather than in the initial state. We then get a Markov chain on the probability 
space of the two points at the ends of the path, 

n x x fly, (22) 

the transition probability depends on the state fj,, but otherwise obeys the properties of a bistochas- 
tic map. This allows the continuum limit of (|20| ) to be taken without leaving elementary probability 
theory. In this limit, we define the densities 



p(x):=ma 3 N X ] vj := a 3 U X , 



and get a damped exponential 



TT (1 - Nx + s'ae) ~ exp { - — / p(x + r'e)dr' 1 , 



0. 



We shall use an identity similar to (pi]) to find the mean free time, subject to two refinements. 
First we require that the site x + sak be empty at the time t(s') that this point is reached by our 
travelling particle. By the same argument, the free path is r = sa if the site x + a(s + l)e is 
occupied at the time at which the arriving particle reaches it. This refinement leads to a transition 
rate is non-local in the time; it will turn out that in the model we construct, we can replace these 
intermediate times by the current time with error of second order. With this done, the transition 
probability is invariant under time-reversal. It then follows from Lemma 1 that, by adjusting the 
stay-put probability so that the rows add up to one, we get a bistochastic map on the two-point 
sample space. 

The second refinement comes from the requirement of (/-invariance in the continuum limit. 
We claim that the 'thermalised part' of a state should be a (/-invariant concept, and this will be 
achieved by (|45|). The Euler dynamics, which is close to the true dynamics, brings the state Jl 
out of LTE in any time interval. The part of the state not thermalised soon becomes thermalised 
by collisions, and it is this thermalisation that is involved in the dissipative part of the dynamics. 
There is some ambiguity in the choice of splitting, because all the particles leaving x seem to be 



instantly dethermalised unless \x is at equilibrium. However, inasmuch as the gradients are small, 
some part of the thermal state at time to might remain thermalised at time dto later. How do we 
decide on how much? It is easy to agree about the rate at which particles thermalise at x at time 
to- We include those particles arriving at x having a momentum k and a relaxation time t' , and 
originating at x — kt'/m. The rate at which they transfer mass, momentum and energy is 1/t' 
times the mass, momentum and energy they carry. 

We note that the concept of thermalisation at a point x over a time-interval (0, t) is not a 
(/-invariant concept. For suppose that in the inertial frame O, various particles thermalise at x 
at times < t\ < . . . < t n < t; then in an inertial frame O', moving relative to O with velocity 
V, they will thermalise (at the same times) at the points x + Vti. We therefore must consider 
the thermalisation (and dethermalisation) that occurs at t in a time-interval (t,t + dt), for an 
infinitesimal dt; this has an invariant meaning. The number thermalising in any time interval of 
length dt is dt times the rate at which the thermalisation occurs. 



The continuum analogue of (20) is a function w(x, k, to; t) which is the probability density that 
a particle at x with momentum k will travel exactly a distance r = \k\t/m and then thermalise at 
y = x + kt/m in the time interval (to + 1, to + 1 + dt). Let W denote the probability density that a 
particle at x with momentum k has had no collision up to the point y, and let C(y, k,t + to) (for 
collisions) be the probability density that a particle at y with momentum k at time to + 1 will be 
thermalised in the tube of diameter a and length dr = \k\dt/m. Then we have 

w(x, k, t ; t) = W(x, k, t ;t)C(x + kt/m, t + t). (23) 

Then the analogue of (|2l|) is 



poo 

/ w(x,k,t ;t)dt = 1 (24) 
Jo 

for all x, k and to- The mean relaxation time is defined to be 

/ w(x,k,to;t)tdt = ti(x,k,to). (25) 
Jo 



We note that <M) 

— POO 

/ W{x,k,t ;t)C{x + kt/m,k,t + t)dt = 1 (26) 
Jo 



can be solved in terms of C: 



W(x,k,t ;t) = exp|- J C(x + kt 1 /m,k,t + t 1 )dt 1 ^ . (27) 
To see (p7|), differentiate to get 



; ;;; - 

and this is just minus the integrand in (|26|). Then one verifies (p(i|): 

/ W(x,k,t ;t)C(x,k,t + t) = - dt = -{W(oo)-W(0)) = l 

Jo Jo ot 

assuming that p is bounded away from zero along the line x + kt/m; this is enough to ensure that 
the mean free path is finite. By construction, w = —dtW and another form for the mean free time 
is 

/■oo 

ti(x,k,t )= W(x,k,t ;t)dt. (28) 



o 



3.3 The Collision Function 



We now find the collision function C appropriate for a hard-core gas. In the dynamics of the lattice 
model, particles hop from site to site with various rates. To implement invariance under Q , we 
extend the fields p, E, vj from A to its convex hull R 3 as continuous piecewise linear functions, 
and consider a particle at x with momentum k not necessarily lying along a basis vector. Although 
the path of such a particle might not intersect A, we must assign a rotated version of w as the 
probability of a free time t, such that ( |24| ) holds. We can assume that the particle hops to the 
nearest site of A to x+kt/m. The continuum limit then makes sense, with densities replacing 
probabilities. Now divide the event, 'the free time is V into the subevents 'the free time is t, and 
A collides with a particle B of momentum q\ Let w(x, k, q,t) be the probability density for this. 
Then 

w(x,k,t) = J d 3 qw(x,k,q,t). (29) 

If t increases to t + dt, then the free path r increases to r + dr, where dt = mdr/\k\. The probability 
that A meets B must be independent of the Galilean frame of reference. Consider the frame O' in 
which B is at rest. Let 

(x,t) i — ^ (x' ,t') = (x + qt/m,t) 

k k - q (30) 

be the Galilean transformation, and denote by p', u' , Q',p' the C- A^-J" variables and the probability 
as observed in O' . Then 

p'(x',t') = p(x,t) (31) 
Q'(x',t') = e(x,t) (32) 
u'(x',t') = u(x,t) — q/m (33) 
p'(y',q',t') = p(y,q,t). (34) 

Then A has momentum k — q and in time dt (which is the same in all Galilean frames) A sweeps 
out a region of volume dV = a\k — q\dt/m, where a is the cross-section. It meets a particle in this 
volume having momentum with probability 

m- l dVp'(y',t)p'(y',0,t). 

By invariance, this is also the probability of collision in the original frame, which is therefore 

a\k - q\(dt/m 2 )p'(y',t')p'(y',0,t') = a\k - q\(dt/m 2 )p(y,t)p(y,q,t), (35) 

by ( |3~l| ) and ([34|). This suggests the choice of collision term 

C(y, k,t + t) = -^ [ d 3 q\k - q\p(y, t + t )p(y, q, t + t). (36) 

Note that we include collisions in which the particle A is hit from behind by the particle B. 

3.4 The mean free time 

Our formula (^6[) for C gives for the mean free time 

tAx,k) = / dtex.pl tt I dt\p(x + kti/m,ti) 

Jo I m z Jo 

J d?q\k - q\p(x + kh/m, q, . (37) 



Comparing this with the the identity 

ran 

rite - */*' 



we see that if the integrand in the exponential in ( J37| ) had been independent of t±, then we could 
have identified ti as 



*t = v {p( x ) J d3 i\ k ~ v\p( x > <?)} • ( 38 ) 



Since e ~ 10 , only values of t less than 16i^ contribute significantly to the integral J ...dt in 
(37). By the mean- value theorem, we may write the argument of the exponential in ( |37| ) as —t/t^, 
where ti is the expression ( |38| ) evaluated at x + kt^/m for the intermediate value £3 < 16i^. Then 
(given that = 0{tq)) the correction to (|38|) is 0(tj), and so (|38|) can be taken as the mean free 
time. Now, p = p + 0(te), so we may put p = p in (|38|). Let 



= y d 3 q\k - q\p(x, q)/p(x, k) 



In the ratio R, the partition function cancels. The exponent in the Maxwellians is 

-TT-q-q-q-C + ^-k-k + k-c = 

2m 2m 

- A(q - fc) • (q - fc) - (g-fc).((/tym)fc + C). 

We change the variables of integration to q — fc, rewritten q, to get 

R = / cPg q exp <^ q • — k + C 

J I 2m \ m 



Put 



thus, ck is the peculiar velocity [14], p. 27. Let q' = (/?/m) 1 / 2 q. Then 

|q|d 3 g = (m/(3) 2 \q'\d 3 q'. 
Dropping the dash, and choosing the axis along k, we get 



R = J d 3 g|q|!ZLexp j-^-g 2 - q . K | 

= / q 3 dq / sin9d6 / dy?— ^- exp < — -g 2 — gK 
j io Jo P 12 

= 2tt / q 3 dq-- r e-< 1 /2 {qn)- 1 (e qK - e~ qK ) . 
Jo p z 



COS ( 



Thus, 



R= 2 -^e^ 2 (I 2 (- K )-I 2 ( K )) 



where we use the functions I n fl7|l 



i n ( K ) = I e -i/2(«/+« 



Put 

F(k) = K exp{--U 2 } (J 2 (-«) - hi*))' 1 . 



2 

Then the mean free time is given by 



ppte = ^F(k). (40) 



3.5 Galilean invar iance 



We now show that the total transition probability rate is invariant under Q. We have split up the 
time evolution into sub-processes, in one of which particle A at x with momentum k has a free 
path of length r = sa, and then collides with particle B of momentum q between r and r + dr. The 
rate of this process was taken to be the same as that of a process in which B was brought to rest 
by a change of inertial frame, and makes a collision between r' and r' + dr', the positions as viewed 
by the observer moving with B. This is obviously necessary if the theory is to be Q- invariant. We 
now show that it is also sufficient: the rate of the one physical process, as viewed in two relatively 
moving frames, will be shown to be the same. The key is to remark that the time interval dt in 
which A collides with B after its free path is the same in all inertial frames, in contrast to the 
distance gone, the free path r and its increment dr, which depend on the speed of A. 

Suppose that O, O' are two inertial observers, with O' moving with velocity — v relative to O, 
such that t' = t and 

x' = x + vt. (41) 
The field variables as viewed by O' are p', p' = ma~ s N', k' , q', where 

p'(x',k',t) = p(x,k,t) 

N'(x',k',t) = N(x,k,t), (42) 
but along a path of a moving particle, the probability densities w and w' must satisfy 

w(x, k, q, r)dr = w'(x', k', q', r')dr'. (43) 

Here, 

k' = k — mv q' = q — mv. 

The point x + kr\ = x + kt\jm on the free path, is assigned the coordinate x+k't\/m by O' . The 
righthand side of (^) is calculated by C using (36) and (|27| ) to be 

eX'+k r' J r l 

J\ 1 



-vfr x a p / (a ; / (^ 1 ),^ 1 y( ^c , (^ 1 ),q',^ 1 )|fc , -g'ljjTj 

N'(x', 0)p'(x', k', 0W (x' + feV, t)) \k' -q'\p' (x'(t), q' , t) , , 

^ / m\k 



dr' 



Here, r[ = \k'\ti/m. Then, by using ([4^) and the remark that 

dr[ dt\ dn 

TT7T = = TTT J < ri < r, 



we see that O' and O assign the same probability to every event, so the integrals over k, q are also 
the same. 



4 Compressible Navier- Stokes with Temperature 



4.1 The Fundamental Relation 

The number of particles, thermalised plus unthermalised, is conserved locally; that is, in any region, 
however small, the loss in particles is the same as the integral of the current over the boundary. 
This local conservation law does not apply to the thermalised subset of particles. A particle can 
cease to be thermalised at x and collide at y, thereby returning to the fold after a time in the 
unthermalised state. We now show how the total probability fx is related to the thermalised part, 
~p,. At time to, any particle at x of momentum k must have been from a thermalised sample at 
some earlier time, to — t, at the point x — kt/m, and remained unthermalised at x, which it passes 
at time to- The probability of thermalising exactly at x is zero. It must thermalise at some later 
time, say after it has travelled for a free time r = t' . Then t' > t must hold. We first compute the 
probability arising from a hop of fixed size r' = \k\t' /m. The rate at which this occurs is 

(l/t')w(x - kt/m, k, to - t, t'). 

In the interval of time from t to t + dt, the number of hops is rate x dt, so the probability of a 
particle being at x at time to with momentum k, and having free time t', is, at time to, 



P(t') := PTob fi {Lo:jV x (u) = lnV x = k\T = t'} = 
rt' _ 

= / dt{t'Y l N{x - kt/m,t -t)p(x - kt/m,k,t -t) 
Jo 

w(x — kt/m, k, to — t, t'). 

Now, w is a density of probability (of collision) as a function of free path size t' , so the total 
contribution to \i due to hops from one side of x to the other along the line of k is 

/■oo 

N(x,t )p{x,k,to) = / P(t')dt' 

Jo 

f°° dt' /"*' 

= / — / dt N(x — kt/m, to — t)p(x — kt/m, k, to — t) 
Jo v Jo 

w(x — kt/m, k, t -t,t'). (44) 

This is the fundamental relation. The logarithmic divergence at t' = is only apparent, if the 
functions entering the integral are smooth enough. If so we can expand in Taylor series in t' up to 
0(ti) around the point x, k,to, at which the functions are evaluated: 

Np = jf* dt \Npw(t') - t(k • d/m + do) (Npw{t')) } 

= Np-\ (— + d )Npt e (45) 



2 \ m 



because of (24) and (|25|). Here, do means d/dto- This is a ^-invariant splitting; for, the equation 
shows that Np and Np differ by a quantity of order ti in smallness, so Np transforms correctly 
up to first order. But ti is (/-invariant, and k • d/m + do is a (/-invariant operator (on fields that 
transform correctly), so Np transforms correctly up to 0(tf); and so on. 

The second term in ([l5|) is responsible for the dissipation. Putting in (^) for tp, we see that 
the density cancels; so the conductivity and viscosity of a gas are independent of the density. This 
is Maxwell's famous result. 

Taking the expectation values of x = "P or £ gives us the relation between the mean of the 
thermalised particles and the true means of all the particles. In this, we need to evaluate do applied 



to the thermalised variables. Here we can assume that these obey Euler's equations, since these 
hold up to first order, and the operator do acts only on small quantities. The integral over w in 
(44) acts as a smoothing operator, so we expect \i to be more regular than /Z. 



4.2 The Euler equations 

The current of the conserved variable x is 



j x ■= J d 3 kNpTx, where T = V/m. (46) 



This gives us the dynamics 

d (x) _ a f ,3., ,fu\ATf„\„f„ U\ k 3 _ A U, k i 



dt 



-d 3 [ d i k X (k)N(x)p(x, fcpl = -S,-(x% (47) 
J m m 



These equations, for x running over the slow variables m, k and k • k/2, can replace C-N-S-T for 
dilute inert gases. Putting x = m gives the usual 'equation of continuity' 

^ + dj (u jP ) = 0, (48) 

which is exact. We shall show how to compute the other equations up to order tg for our choice of 
C. Our strategy is to us e (|45| ) in (f47|), allowing us to evaluate the righthand side in terms of the 
means in ~p; we then use (|4q ) again to rewrite this in terms of the true means. 

The zero th order approximation to (|45|), namely Np = Np, can be put in (|47]) and computed 
exactly: we get the Euler equations; this is shown very smoothly by using the cumulant generating 
function, log Z. 

For the momentum, put x = T^x- Then we have for each x, 

d{Nut) - - dj (N^{TT^ 



dt 

= -dj (N({T T j ) T + «*« 



Since 

we get the Euler equation for momentum conservation: 





dt (pu l ) + djipu'u 1 ) + di(pk B Q/m) = 0. (49) 



Finally, for the energy, put x = &x = k • k/(2m) in (|47|), which then becomes 

E = -div (m^NEplSV]) . (50) 

Now, for each x, 

Bp[£V j ] = (£r j ) T + (£)(v j ) 

d 2 log Z 



d(3d(j 



+ NEmv? 



Tn^ ' ■ f 1 
— + mNu 3 [ k B Q H — u • u 



p 2 2 



from (||). Since 

= -(3u j and E = mN(3k B 8/m + u ■ u)/2 
we can collect terms to get the Euler equation for energy conservation: 
d 

— {p(3k B e/m + u-u))/2 + div (pu(5k B Q/m + u-u)/2) = 0. (51) 

The pressure appearing in the usual form of Euler's equations is here replaced by pk B @/m, the 
pressure for a perfect gas. This differs from the static pressure (|l7| ) by terms which vanish in the 
low density limit. We shall use the Euler equations, which are first-order PDE in space and time, 
to relate 8q to a first-order gradient. Let D := u • d + do be the Lagrange material derivative. 
Then we have 

Lemma 2 The Euler equations §Mi), (4&) an d ( PjU imply the short Eulers: 

Dp = -pdju j (52) 
Eh? = -k B p~ l di{p@/m) (53) 

DQ = --Gd jU j . (54) 



Proof. 'It does not seem necessary to reproduce the details of this proof; the mathematician will 
be able to construct them for himself, while the physicist will probably not wish to be detained 
over them' 

The <5-invariance of this form is obvious. 
4.3 Calculations 

It is the full state p, rather than the thermalised part, /Z, that is usually measured in experiments. 
For example, a measurement of density can be made by noting the absorption of a laser passing 
through the gas. The scattering of photons with particles makes no distinction between thermalised 
and non-thermalised particles. Again, one can measure the temperature by a probing thermometer, 
which would tend to thermalise any particles that struck it, whether they were thermalised before 
or not. However, the means in the thermalised state are much easier to calculate; the state p is 
Gaussian, and the fields are independent at different points (at the same time, say to)- This enables 
us to relate the extensive to the intensive variables. For fj,, the intensive variables have not even 
been defined yet. 

In the Boltzmann equation, authors write the phase-space density as a product N(x)f(x, k) 
with some hesitation, as and "P are not independent random variables, even in the Maxwellian. 
Not to worry. For the general state p we define 

N(x) = Prob M {uj : N{x) = 1} , (55) 

and p(x, k) is the conditional probability 

p(x,k) = Pmb^{V x = k\M{x) = 1}. (56) 

By Bayes's definition, 

Px(ljx) = N(x)p(x,k). (57) 

Then we may define p = mN/a? . The momentum density also has a definition in terms of p, which 
does not assume that p is Maxwellian, by 



zv = Il/a 3 = E^[7']/a 3 . 



We can then define the velocity field, without recourse to information geometry, by 

u = zu/p, 

provided that the density never vanishes. This is indeed so, as we see from the basic equation (fl4|). 
Finally, we define the thermal energy per unit mass, e, to be given in terms of a 'temperature' G 
at each x by 

3 1 
e = —k B Q := m^BplSx] - -u ■ u. (58) 

This definition of temperature is (/-invariant, as it can also be written as mc 2 E p [KjKj)/2]. 

Suppose that we know the fields p, and it, at time=tcb and thus also their space gradients. We 
can use (45) and the short Eulers to find the barred parameters of the Maxwellian p, to first order 
in t£. When any expression is multiplied by tg, we are able to replace any thermalised parameters 
by the above unbarred parameters, or vice versa, with only a second order error. So the analysis 
reduces to linear algebra. 

For any local random variable x, slow or not, a~ 3 x is its density, and 

a~ 3 (x) := J d 3 km~ 1 p(x)p(x,k)x(k) a~ 3 (x) := J d 3 km~ l p(x)p(x,k)x(k). 

To help in the evaluation of various derivatives of integrals that occur here, we note that while F(n) 
is a complicated function of u and (5, the integrals arising can be evaluated if we change the variable 
of integration from k to k, given in (|3^). We do this for each x, and it is valid provided that we 
keep the derivatives d and do to the left of the expression. We note that d 3 k = {mj ji) 3 ! 2 d 3 K. Then 
using ( [45| ) and ( |40| ) we have to order if. 



a-'\(x) ~ <X» = I d 3 k (— + do) F( K )f3 2 X (k) 

1 / m \ 1 / 2 f 

— \^—J dj@~ l/2 J d 3 K(cKj + Uj)x[m(cK + u)]F(k) 



, 1/2 

AlTCT 

1 / m\ 1 / 2 



4"7T<7 \k B 



(jp) d @- 1/2 J d 3 Kx[m( C K + u)]F( K ). (59) 



Let us put 



m I m 



D 



1/2 



\ n ■- -Jl I £L ) / K 2n F(n)dK, n = 1,2,3. (60) 



On putting x = m, the term odd in k is zero, and we get 

p = p- X 1 d j (®- 1/2 u^ - Ai^oG- 1 / 2 . (61) 
We now put x and = (sin cos 99, sin sin 93, cos 93). In the following calculations, we use that 

/ dQkikj = I sm9 d6 dip kikj = (62) 

and that the integral of odd powers of k are zero. Then ( |59|) gives 

1 k £ 
3 m 

We now put x = & * k/(2m); we get 

1/2 



Wi-Wi = -\^X 2 d l Q 1 ' 2 - X 1 d j (e-^UiUj) - Xido (@- 1/2 Ul ) . (63) 



a~ 3 (E -E) = — <9 O 0~ 1/2 / F(n)K 2 dKdn (at* + uj) (ckj + uA 

8ttu \k R J J 

djQ^ 1 / 2 J n 2 dK,dfl F(k) (cKj + Uj) (cKg, + ug) (ck^ + ug) 



in ! m ^ 1 ^ 2 



87r<r \k B 



which simplifies to 



k B X 2 „ „i 10 A 



d e 1/2 - Y d ° (@~ 1/2 w) . (64) 



2m 2 

We can use these results to relate (x) and (x) for any polynomial. Let 

$X = (x) - (X>- (65) 
Then up to 0(te), <5 is a derivation, and we have for example, 

5(pu • u) = 2uj5(puj) — u • u5p. (66) 



Lemma 3 We have the relation 



5(pQ) = -\\ 2 Q l ' 2 d jUj - \m\2D® 1 ' 2 = -^Q^djUj. (67) 
y o i) 



Proof. 



h B 5( P Q/m) = |^(pG-p6) 
2 2 m 



a 5E — Uibwi + UiUi5p/2 
6m 2m 



making use of (pip, (S3) and (p3). This gives the result, using Lemma (0) to get the second form. 



4.4 Equations of Motion 

The equations of motion express the rate of change of the slow variables (the conserved quantities) 
in terms of the mean of their microscopic currents; we use 



4.4.1 Viscosity Terms 

Let Jij. be the j' th component of the mean BBGKY current density of the z th component of mo- 
mentum in the state p, and let J^. be the Euler term puiUj + k B 5ijpQ/m. Let J 3 m . be the same 
object averaged in the state p; as found in (]49|), this is the Euler term with the barred values of 
the fields. Then the choice x = hi in (^7j) gives, using (|59[), 

m? 

7 W . - £^ (j^J 1 8 e- 1/2 J K 2 F(K)dKdn(cKi + Ui)(cKj + 



J4 = / afk^pp 



, 1 1I 

m ( m \ >~ 



47TO" V k B 



8e@ 1 ^ 2 / K 2 F(K)dK (cKi + Ui){cKj + Uj)(cK£ + Ul). 



Noting that the odd powers of Kj integrate to zero, this reduces to 



JL 



y 

*J -7-TT 



k B \2 
3m 



a t (e 1 / V) 8ij + dj (e 1 / 2 ^) + di (e 1 / 2 ^) 



UiUjUg J + 



V2, 



dm 



(68) 



Note that derivatives act on all functions to their right. We get our momentum equation by relating 



JL to rj,-. 



k r 



J m = J% t ~ Uj5(wi) - Ui5(wj) + (uiUj5p) -8(pQ)5ij, 



in 



(69) 



for which we use (61), ( )63| ) and Lemma [3|. Putting this in (|68|), we see that most terms cancel, 
leaving 

j£. = J°J t - \^M^ 1/2 Uui + diUj - \d t u E Sij) . (70) 



The kinetic pressure is traceless, so there is no bulk viscosity [35]; we get a very special case of the 
general equations, at the edge of possible values. It could be that the Stokes relation only occurs 
when, as here, there is no interaction whatsoever between the particles outside the hard core. It is 
interesting that Ai does not occur in the answer, and that all the terms in do cancel out without 
recourse to the short Euler equations, Lemma (|2[). The viscosity coefficient increases as G 1 / 2 with 
temperature, like Enskog's and Chapman's prediction from Boltzmann's equation. 



4.4.2 Equation of motion for the energy 

Let J 3 be the mean energy current density, with a bar if the state Ti is used, and let 



5k Br, 1 

fc> H UiU; 



,2 m 2 

be the Euler energy current. Putting x = kiki/(2m) in Q47|) gives 



where 



iv 



e-i/2 



8na 



J 3 =f E+ d [A]+d e [B], 



m \ i/ 2 f 

— J J K 2 d,K dttF(K) (cKi + Ui) (cKi + Ui) (cKj + Uj 



^l n -l/2 5/c s A 2 i/2 

B ' UiUiUj 6 ' Uj, 
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(71) 
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(73) 
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k » x % e e*/ 2 - 7 J^le^ UjUe - ^5^ 2 u,u, - ^e- 1 ^ 

6m 
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6m 



UgUiUiUj. 



(74) 



Now, J E — J° 3 = —8J 3 E has six terms: 



5k B 
2m 



(@5wj + Uj5(pQ) — @Uj5p) 



+ 



—UiUidvuj + UiUjdvui — UiUiUjbp 



(75) 



This can be evaluated using (|6l|), ( |63j) and Lemma (|3|). Collecting up A, B, and —SJ^, using the 
same method of proof as in Lemma (|2|), the heat equation reduces to 



,dE 



K fX 



5A9 5Ai 



dt 3 E m V 4 4 2 ' J J 



e 1/2 Ui d iUj 



@ 3/2 djlogp 



The coefficient of the Fourier term is positive, since 

1 



2k r 



X 2 dj 



® l l 2 UjdiUi 



e 1/2 d jU iUi/2 



(76) 



+ — k > 0. 



The new term involves the logarithmic derivative of the density, whose presence in a gas of a single 
component is denied in the literature. It means that a gradient in the density contributes to the 
heat current. This is the Dufour effect, also called the diffusive thermal effect. The sign of the 
term does not need to be definite. Some authors invoke Onsager symmetry to eliminate this term 
without having to evaluate it, since its Onsager dual, the Soret effect in the continuity equation 
for the mass, is absent. However Onsager duality is not true here, because the state p is not in 
LTE and the transition rate depends on x. The present work does suggest that the effect should 
be looked for experimentally, in say Helium, but this is quite delicate since the Dufour effect is 
transitory, and becomes less pronounced, and is masked by heat conduction and convection, as time 



goes by. The other terms have appeared in the literature [15]. One can check that the system of 
equations is invariant under Q: 'in the tradition of British applied mathematics, it is not considered 
gentlemanly to press a colleague for a proof (G. Pistone). 



5 Conclusions 



We have shown that the 'method of Maxwell' [29] can be made (/-invariant, and gives C-N-S-T 
with a Dufour term. The fluid equations we get are the following: 



d_ 

dt 



dp 

~d~t 

dwi 

~~dT 

[p(e + UiUi/2)} 



-dj(puj) 

-dj(puiUj) — diP + Xdj 



o 1 / 2 f ana + a 



-dene 5, 



-dj ( pUj (e + um^) + ujP) + x 4 dj (e^fye) + x 5 dj (e^dj i og( 



+ xdj 



Q l l 2 Ui (djUi + diUj — -d£U£ 5{ 



(77) 
(78) 

) 

(79) 



The transport coefficients are independent of density, as found by Maxwell. Our starting point is 
not the Boltzmann equation, but a non-local integral equation, (f44|). The presence of coefficient A5 
is at variance with the results of Chapman and Cowling [Q . The Boltzmann equation suffers from 
the Hilbert paradox ||], II, p. 348. Namely, the state in Boltzmann's equation is parametrized 
by the initial distribution /, which is a general integrable function of six variables, whereas the 
hydrodynamic solutions are parametrized by five functions of three variables; the set of hydrody- 
namic solutions cannot describe the most general solution. To show that they nevertheless provide 
a good approximate solution is an extra burden if one uses the Boltzmann equation as the starting 
point. In [14] this question is discussed but not solved on p 120; the authors refer the reader to 



p8|| . Chapman himself has said that reading his book is like 'chewing glass' fl3|| . The fact is that 
the Boltzmann distribution / is too detailed a description for an easy move to thermodynamical 
variables. Hilbert's paradox also shows up as follows. If the initial state happens to be in LTE, 
then the collision term in BE is zero, and (at that instant) the thermodynamic variables follow the 
Euler equation, and the instantaneous rate of entropy production is zero. However, fluids following 
the equations supposedly derived from the BE do not at any instant follow the Euler equations, or 
possess a zero rate of entropy production, except in the special states with © and u independent 
of x. This paradox lasts a very brief time, after which the stirring due to the Euler convection 
spoils LTE. For consistency between the BE and the fluid equations, we are not allowed to choose 
an arbitrary initial state for the BE; the small deviations from LTE must be related to the fluxes 
of the theory Q, p. 160. Our approach avoids the paradox: an initial state fi can be LTE, and 
then we would modify (^4|) in the obvious way. 

The point of view of the present paper differs from the usual one, such as [[^J. We assume 
that just after a collision, which we prefer to call a thermalisation, the particle is well described 
as being in LTE, and almost independent of its neighbours. This is a good time 'to carve Nature 
at it joints' |2j, p 341. It is during the free propagation that the state loses its LTE property, 
since then particles from regions of different density and temperature come together. They are 
likely to be independent because 'they have a different history' ||34f| . Concerning the other phase 
of the dynamics, our collision term involves the density at different space-time points, whereas the 
Boltzmann kernel is local. In the models in |3^| , several collisions are needed before a particle is 
close to thermalised; the authors do, however, show a preference for models, called 'efficient', in 
which the collision output is randomised over the available channels, and so thermalises rapidly. 
In the present model, a particle thermalised after one collision. Our result shows that the details 
of how many collisions are needed does not affect the qualitative results, though it changes the 
relationship between geometric cross-section and mean free path. It is difficult to believe that the 
simplification made here is responsible for the Dufour effect. 

Our equations ( |79| ) may be less singular than C-N-S- T since the presence of the Dufour term 
means that the symmetrised part of the operator 8q has a principal symbol of full rank, at least in 
general position. It might be that C-N-S-Thas no global smooth solutions, or, even if it has some, 
it might be too hard for anyone to prove it. Whatever the case, it might be easier to show that 



(44) has smooth global solutions 

We can generalise in various ways. If there is an external potential, it does not affect the 
local state, because it cancels out in p,; however, it does affect the hopping rates, and thus appears 



in the equations of motion. In a paper [19| we find the equations of motion for a fluid moving in 
a potential, in a non-galilean model. The same method can be applied to the present model. It is 
possible to extend both models to the case of inter-particle potentials by following a suggestion of 
Biler and collaborators J?], ||, |9|, |l^, 11, 12, 30 1. This gives a macroscopic dynamics in which the 
rate of change of energy at a point x is governed by the mean field of all the other particles. It 
seems unlikely that making the model more realistic by including interactions will exactly cancel 
out the prediction of the Dufour effect; thus this should be looked for in Helium or Argon. 
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